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Abstract 

The properties of hot matter are studied in the frame of the rel- 
ativistic Brueckner-Hartree-Fock theory. The equations are solved 
self-consistently in the full Dirac space. For the interaction we used 
the potentials given by Brockmann and Machleidt. The obtained 
critical temperatures are smaller than in most of the nonrelativis- 
tic investigations. We also calculated the thermodynamic properties 
of hot matter in the relativistic Hartree-Fock approximation, where 
the force parameters were adjusted to the outcome of the relativistic 
Brueckner-Hartree-Fock calculations at zero temperature. Here, one 
obtains higher critical temperatures, which are comparable with other 
relativistic calculations in the Hartree scheme. 
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I Introduction 



The properties of hot and dense nuclear matter play an essential role in the 
understanding of high-energy heavy-ion collisions, supernova explosions and 
proto-neutron stars. For that reason the problem of hot nuclear matter has 
been studied over the last decades in several investigations, which however 
were predominantly performed within the nonrelativistic scheme using 
either effective density dependent interactions 0, H, 0, Q or the Brueckner 
approach [|I], |^, ^ |^, |]. In the relativistic approach investigations of the 
equations of state for T 7^ are relatively scarce. The majority of such 
calculations were performed in the relativistic Hartree approximation (RH), 
where the extension to finite temperatures is straightforward. Details of 
this scheme are given, for instance, in Refs. 0], [0|-|T^. More complicated 
are the relativistic Hartree-Fock- Hl^ ITSU and the Brueckner-Hartree-Fock 
approximation 



T6| , and the application to finite nuclei [0. 



In this contribution we will concentrate on the relativistic Brueckner- 
Hartree-Fock treatment (RBHF) of symmetric and asymmetric nuclear mat 
ter generalizing the formalism as described in Refs. [[ 
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to T 7^ 0. The 

RBHF-approach seems to be of special interest, since it is known for T = 
that the resulting EOSs are much stiffer than their nonrelativistic counter- 
parts . To our knowledge such an investigation has been only performed so 
far by the Groningen group for symmetric nuclear matter |jl6| . As described 
in more details in Refs. UI 



their method uses a nonunique ansatz for 
the T-matrix in terms of five Fermi invariants, which can lead to ambiguous 

1). In order to avoid this 



results for the self-energies (see, e.g., Refs. ||2T 



problem we solve, according to the original scheme of the Brooklyn group 
||23| , the RBHF-approximation in the full Dirac space, which is more tedious 
(for a more detailed discussion, see Ref. [@). Since the formal structure of 
the problem is the same as for T = 0, where one has to solve three coupled 
equations, namely the Dyson equation for the one-body Green's function G, 
the (reduced) Bethe-Salpeter equation for the effective scattering matrix T 
in matter and the equation for the self-energy E, we will not repeat here 
the equations. As in Refs. ITB, IT9| we will restrict ourselves to the incor- 



poration of intermediate positive-energy nucleon states only, where now the 
Fermi step functions are replaced by the Fermi distribution functions n^{T) 
at finite temperature T (for details, see Ref. p5|). The Green's function obeys 
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for T 7^ the spectral representation 

J [Pp-UJ-lT] Pp-UJ + lTTj} 

with 

= + ^ e(-a;) . (1.2) 

The formal structure of the spectral function A{'p,uj) [|TB| remains unaltered 
to the case for T = 0. 



A further difference in comparison with Ref. |]T6[ is that we take the mo- 
mentum dependence of the self-energies into account. Since the pole of the 
quasi-particle propagators occurs for T 7^ in the integration domain of the 
intermediate states, one obtains, in principle, complex effective scattering 



matrix elements and self-energies. It was checked in Ref. [^] that the imagi- 
nary part of S turned out to be small. Therefore we neglect also Im S in the 
calculations. For the one-boson-exchange interaction we used the modern 



potentials constructed by Brockmann and Machleidt [^. We select for the 
presentation the so-called potential B, which gives the best results for the 
EOS {E/A = -15.73 MeV; po = 0.172 fm'^; K = 249 MeV; J = 32.8 MeV) 



at zero temperature in RBHF-calculations (see Refs. |T^, [T9[ , for the poten- 
tial A the outcome is similar, see Ref. [§^]). For the sake of comparison we 
also treated the RHF-approximation, where we adjusted the force parameters 
to the outcome for the EOS for symmetric and asymmetric nuclear matter 
at T = E^. The RHF-approximation has in comparison with the 



RH-approximation the advantage to resemble in its formal structure more 
to the RBHF-approximation with the benefit of a much simpler numerical 
treatment than in the RBHF case. 

For finite temperatures one needs for the determination of the pressure 
the free energy per baryon, defined as 

f = u-Ts, (1.3) 

with the entropy per baryon: 

Y, j <CP Hp) In n(p-) + (1 - „(p-) 1„ (1 - nm ■ (1.4) 
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The pressure is given by: 



(1.5) 



T,p-r 
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II Results and discussion 



Common to almost all nomelativistic treatments are EOSs typical in shape to 
the standard Van de Waals behavior. The value for the critical temperature 
depends strongly on the choice of the forces (and approximations). The 
bandwidth reaches from 14 to 22 MeV, with a critical density of about 1/3 
of the saturation density po [H • In a recent calculation of asymmetric matter 
for T 7^ a critical temperature of 20.8 MeV (pc = 0.39 po) was obtained for 
symmetric matter, which decreased to 8.0 MeV (pc = 0.24 po) for the limiting 
asymmetry of 5 = 0.75 0]. In the relativistic treatment within the framework 
of the relativistic Hartree approximation the Van der Waals behavior is still 
present, but in general one obtains lower critical temperatures (for instance, 
T, = 14.4 MeV, p, = 0.31po [0, T, ^ 14 MeV 

Of great theoretical interest is the microscopic RBHF-treatment for two 
reasons: Firstly, as mentioned before, the RBHF-EOSs are much stiffer than 
their nonrelativistic counterparts ||2^. Secondly, according to Refs. |jl], |16 



the typical Van der Waals behavior may be questionable, and one obtains 
a rather low critical temperature, Tc ~ 8 — 9 MeV |I| or 12 MeV (pc — 
0.6 Po) [|16|) depending on the approximative treatment of the self-energy 
(see, also the discussion in Ref. [0). Due to these reasons we recalculated the 
RBHF-approximation in the full Dirac space as described in Refs. [13, 123 



with the Brockmann-Machleidt potentials which give better results for 



the saturation properties [18|. As discussed before we treated also the RHF- 



approximation, where the force parameters are adjusted to the outcome of 



the RBHF-treatment at T = 19, 25 



In Figs. 1 and 2 we present first on a larger scale the EOSs for different 
asymmetries and temperatures, computed for RBHF and RHF, respectively. 
On this scale the energy per nucleon differs not very much in both approxima- 
tions. This holds also for the pressure (see Fig. 3), free energy and the entropy 
(see Fig. 4). As expected the EOSs seem to be stiffer as in non- relativistic 
Brueckner calculations (cf., for instance, with Ref. 0). The entropy behav- 
ior is similar to Ref. [0] and agrees with the experimental situation (see, e.g.. 
Fig. 6 of Ref. 0]). Interesting is the EOS in the lower density domain, where 
one expects a phase transition. The situation is depicted for symmetric mat- 
ter in Figs. 5 and 6, from which we deduce a critical temperature of 10.4 MeV 
(RBHF) and 15.2 MeV (RHF), respectively. (For asymmetric matter, see 
Ref. [0.) The RBHF value is in accordance with Ref. |TB| and with one non- 
relativistic BHF - calculation [|^ . The RHF result for Tc is higher and agrees 
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more or less with the result of Ref . \T2\, where one uses a smaller equilibrium 
density. The critical temperatures and pressures for different asymmetries in 
the RBHF-treatment are given in Fig. 7 (for RHF, see Ref. |^). The critical 
densities are shown in Fig. 8. For symmetric matter they are slightly above 
1/3 of the saturation density but smaller than in Ref. |]TB . 



One should however mention in this context that in general BHF-calculations 
are facing numerical convergence problems for small densities and may 
even not be applicable |^ . For that reasons the conclusions drawn from the 
BHF-approximation for low densities (and temperatures) should be consid- 
ered with some caution and a phenomenological description with adjustable 
parameters, for instance, with Skyrme forces, might be a suitable alternative 
in this sensitive domain. Unfortunately the critical temperature can not be 
extracted directly from experiment and additional models and assumptions 
are needed to obtain it from hot nuclei produced in heavy ion collisions 
However more recent experiments imply that a lower critical temperature 
might be possible 



As a final point we would like to mention that the symmetry energy 
remains - as expected - a monotonic increasing function of the density in 
the case for finite temperatures (see Fig. 9), so that the composition of proto- 
neutron stars should show, as in cold neutron stars for relativistic EOSs, the 
tendency to lower asymmetries with increasing density H, |I4|; contrary to 



some nonrelativistic EOSs with a non-monotonic symmetry energy 



III Conclusions 

In conclusion, we have performed a calculation of hot symmetric and asym- 
metric nuclear matter within the relativistic Brueckner-Hartree-Fock scheme 
using modern OBE-interactions constructed by Brockmann and Machleidt. 
It turned out that the critical temperatures are smaller than it is the case 
for the majority of nonrelativistic treatments. We have additionally treated 
the relativistic Hartree-Fock approximation at T 7^ 0, where the Lagrangian 
parameters were adjusted to the outcome of the RBHF-treatment for T = 0. 
Here the critical temperature is in the range of other relativistic treatments 
performed in the Hartree scheme. 
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Figure captions 

Fig. 1: Energy per nucleon for different asymmetries {6 = ^^2— ££) and tem- 
peratures as function of the density in the RBHF-approximation 
(Brockmann -Machleidt potential B). 
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Fig. 2: Energy per nucleon for different asymmetries and temperatures in 

the RHF - approximation. 
Fig. 3: Density dependence of tlie pressure for different asymmetries 

and temperatures in tlie RBHF - approximation. (In tfie RHF- 

approximation tlie pressure increases less for higher temperatures). 
Fig. 4: Entropy per baryon for nuclear matter versus density at differ- 
ent temperatures (for the experimental comparison, see Fig. 6 of 
Ref. 0). 

Fig. 5: Pressure as a function of baryon density for nuclear matter {6 = 0) 

at different temperatures in the RBHF - approximation. 

Fig. 6: Pressure versus baryon density for nuclear matter at different tem- 
peratures in the RHF - approximation. 

Fig. 7: Critical temperature and pressure as function of the asymmetry 

(RBHF; Tc and Pc are smaller than in the RHF). 
Fig. 8: Critical density as function of the asymmetry. 

Fig. 9: Symmetry energy as function of the baryon density at different 

temperatures (RBHF). 
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